Quantum Phase Analysis of Extended Bose-Hubbard Model 
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We have obtained the quantum phase diagram of one dimensional extended Bose-Hubbard model 
using the density-matrix renormalization group and Abelian bosonization methods for different 
commensurabilities. We describe the nature of different quantum phases at the charge degeneracy 
point. We find a direct phase transition from Mott insulating phase to superconducting phase for 
integer band fillings of bosons. We predict explicitly the presence of two kinds of repulsive Luttinger 
liquid phases, apart from the charge density wave and superconducting phases for half-integer band 
fillings. Our study reveals that extended range interactions are necessary to get the correct phase 
boundary of an one-dimensional interacting bosons system. 
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1. INTRODUCTION 



We present the quantum phase diagram (QPD) of ex- 
tended Bose-Hubbard model (EBH) 0,0,(1], for spinless 
bosons on a one-dimensional lattice. These bosons could 
represent Cooper pairs undergoing Josephson tunneling 
between superconducting islands or Helium atoms mov- 
ing on a substrate. Here we are interested in the quantum 
phase transitions of the ground state. This transition is 
driven by quantum fluctuations of the system which are 
controlled by the system parameters [J, 0] • There are a 
few experiments carried out on one- dimensional Joseph- 
son junction systems of different lengths under a mag- 
netic field [!, 0, 0, 0] ; a superconductor-Motts insulator 
transition is observed. In these systems, relevant parti- 
cles are the Cooper pairs which can be treated as spinless 
charge bosons. Hence the Bose-Hubbard model and its 
variants are prominent candidates to study such systems. 
It is also revealed from experiments that the inter-particle 
interactions arc of finite range 0, 0, 0, 0, 0] ■ We study an 
EBH model that incorporates the nearest-neighbor (NN) 
and next-nearest-neighbor (NNN) repulsive interactions 
in addition to the on-site Hubbard repulsion. We also 
consider the NN and NNN hopping of bosons on the lat- 
tice. In what follows, we derive an analytical expression 
for Mott insulating gap at the mean-field limit but this 
alone is not sufficient to capture the effects of all inter- 
actions especially at the charge degeneracy point [§] . We 
map our EBH model to a spin chain model, for strong on- 
site Coulomb repulsions. In this letter, we first present 
the results of Abelian bosonization study of spin chain 
model, we then discuss the results of our Density Matrix 
Renormalization Group method (DMRG) [l(| study of 
the EBH model. 

The model Hamiltonian representing, arrays of super- 
conducting (SC) nano-island, here after refered to as su- 
perconducting quantum dots (SQD) is given by, 
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in the usual notation. 



2. MODEL HAMILTONIANS AND 
CONTINUUM FIELD THEORETICAL STUDY: 



We recast our basic Hamiltonians in the spin language 
IH to obtain; H n = -2 E, n £ 2 (S+S m + h.c), Hj 
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Our total Hamiltonian is, 
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H = Hji + H.J2 + HecQ + HecI + Hec2 



(2) 



Here h — (N — 2n — l)/2 and N tunes the system to 
a degeneracy point by means a gate voltage (eA^ ~ V g ) 
where eN is the average dot charge induced by the gate 
voltage. The phase diagram is periodic in N with pe- 
riod 2. Here we only consider a single slab of the phase 
diagram, < N < 2. Correspondence between the pa- 
rameters of the EBH model and the spin chain model 
is; (n)ti ~ E n , (n)t 2 - Ej 2 , U - E co , Vi - E zl , V 2 
- E z2 12]. When the ratio 0, the SQD array is 

in the insulating state with gap ~ E C o, it costs an energy 
~ Eqq to increase by one, the number of pairs on a given 
dot. Exceptions are the discrete points at N = 2n + 1, 
where a dot with charge 2ne or 2(n + l)e has the same 
energy because gate charge compensates the charges of 
extra Cooper pair in the dot. At this degeneracy point, a 
small amount of Josephson coupling leads the system to 
the SC state. To analyze the phases explicitly near this 
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degeneracy point, we recast the Hamiltonian of the spin 
chain model to a spinless fermion model through Jordan- 
Wigner transformations and carry out a continuum field 
theoretical study of spinless fermion model 1 1311 . Our con- 
tinuum models, in terms of bosonic fields [LJI, are 



Hi = H 
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d x </> dx. (4) 
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Jdx[(d x 9) +(d x (j)) ], is the non- interacting 
part of the Hamiltonian (Eq.2). The velocity, v, of low 
energy excitations is one of the Luttinger liquid (LL) pa- 
rameters while K (Eqs. 3 and 4) is the other. 

<fi field corresponds to the quantum fluctuations 
(bosonic) of spin and 9 is the dual of </>. They are related 
by 4»r — 9 — (f> and <f>L — 9 + (j>. The effect of applied 
gate voltage on the SQD appears as an effective magnetic 
field in the spin representation of these Hamiltonians. In 
Hi, we consider the NN Josephson coupling, on-site and 
NN charging energies of SQD and also the NN charging 
energy correction for the co-tunneling process. In H2, we 
consider all the interactions present in our Hamiltonian 
besides all terms of co-tunneling process [ill , [lit ] . Ana- 
lytic expressions for Ezn and Ez\2 are [Ezi — fjlfe') 
and (Ezi - ff^ - \E Z i) respectively. 



3. NUMERICAL STUDIES 

We use the DMRG method to numerically study the 
QPD of the Hamiltonian in Eq. 1 . We employ the infinite 
DMRG algorithm keeping 128 density matrix eigenvec- 
tors. The site boson Fock space is truncated to three 
and the number of sites in the chain is restricted to 128. 
Per site bosonic filling fractions 0.5 and 1.0 are studied. 
Accuracy of the method is checked by comparing the 
ground state energies, various correlation function and 
charge gap from DMRG studies with exact diagonaliza- 
tion studies of small systems. We have also reproduced 
the result of earlier DMRG calculations satisfactorily Q . 
The discarded density in the DMRG calculations is less 
than 10~ 14 in the charge density wave (CDW) phase at 
p = 0.5 as well as in the Mott-insulating phase at p = 1.0. 
However, the discarded density is slightly less than 10~ 10 
in the SC phase. For calculating the correlation func- 
tion, we target only the ground state. To obtain the 
Berezinskii-Kosterlitz-Thouless (BKT) transition point, 
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FIG. 1: Quantum phase diagram (Eji vs. V g ) of the SQD 
array. A is the Mott insulating phase, Bi is the CDW phase, 
B2 is the Dimer-order density wave, C is first kind of Re- 
pulsive Luttinger liquid (RL1), D is second kind of Repulsive 
Luttinger liquid (RL2) and E is the SC phase. Pi and P2 
are the multi-critical points (please see text). K is one at 
the Luttinger liquid and SC phase boundary and the Mott 
phase and SC phase boundary. K = 1/4 at the phase bound- 
ary between RL1 and CDW state. K = 1/2 and 2 at P2 
(charge degeneracy) and P3 (particle-hole symmetry) points 
respectively. 



we plot the correlation function < 6+ bo > vs r on a log- 
log plot, here r = is the middle site in the chain. The 
slope of this plot in the region between r = 16 and r = 36 
gives —K/2, where K is the Luttinger liquid (LL) param- 
eter (Eqs. 3 and 4). The BKT-point corresponds to slope 
of -1 for p = 0.5 and -0.25 for p = 1.0. 



4. PHYSICAL ANALYSIS OF QPDS. 

Here we analyze the QPD of Hamiltonian Hi (Eq. 3). 
The analytical structure of this model is the same as that 
of the XXZ Heisenberg chain in a magnetic field. Hence, 
LL parameter of Hi can be calculated exactly by using 
the Bethe-ansatz techniques. Analytical expression for 
K [H is of the form 



K 



7r + 2 sin -1 A ' 



(5) 



where A is the XXZ anisotropy parameter of the Heisen- 
berg chain. Our analysis will use the two limits of the 
parameters A, namely Ai and A2; the analytical expres- 
and (2fa - |f£) respectively. 



sions for these are 



We calculate the analytical expression for the LL pa- 
rameter for Hamiltonian H2 (Eq. 4) as 
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In the limit A = Ai, for Eji < 2Ezi and relatively small 
field, the anti-ferromagnetic Ising interaction dominates 
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FIG. 2: 
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Phase diagram of l/£ vs. ii for (a) p = 0.5 and 
1.0 boson density for different values of V2 (£ 2 ~ 



^2 r r 2 < bj.bo >/^2 r < btbo >). We predict three quantum 
phases: CDW, RL2 and SC phase for p — 0.5 but there is no 
evidence of RL2 phase for p — 1.0. Inset in (a) show structure 
factors for p — 0.5 . 
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FIG. 3: Phase diagram in (ti, l/£) plane for boson density 
of 0.5. We predict effect of finite value of t.2 on two phase 
boundaries (CDW-RL2, RL2-SC). 



the physics of anisotropic Heisenberg chain. When the 
field is large, i.e., for large applied gate voltage, the chain 
state is in a ferromagnetic state. In the language of 
interacting bosons, the Neel phase is the commensurate 
CDW phase with period 2. In Fig. [TJ this phase region 
is described by region B\. The ferromagnetic state is 
the Mott insulating state, this is the phase A of our 
quantum phase diagram (Fig. [T]). The emergence of two 
LL phases can be attributed to following reason: RL1 
(K < 1/2) occurs due to commensurate- incommensurate 
transition while RL2 (K > 1/2) occurs due to the 
critically of the XY model. The physical significance of 
RL1 (region C in Fig. [T]) phase is that the coupling term 
is relevant and the applied magnetic field (the applied 
gate voltage), breaks the gapped phase. In the RL2 
phase (region D in Fig. [T]) none of the coupling terms is 
relevant due to large values of K (> 1/2). These phase 
regions, are similar for the two different limits, A = Ai 
and A = A2. The analysis of the physics of two RL 
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FIG. 4: Density-density correlation function for two values 
of V 2 at Vi = 0.4. o - o ti = 0.05, □ - □ ti =0.1 and 
O- O t x = 0.15. 



phases is entirely new. 

In Fig. H we present the DMRG study of EBH 
model (Eq. 1) for boson densities, p, of 0.5 (Fig. [2h.) 
and 1.0 (Fig. [2Jd). We predict three quantum phases, 
CDW, RL2 and SC for p = 0.5. There is no evidences 
of LL phase for integer boson density. So our field 
theoretical analysis is consistent with DMRG studies. 
The structure factor in inset of Fig. [5^ for p = 0.5 
precludes CDW phase for this parameter range. There 
is also no evidence for RL1 phase; one obtains the RL1 
phase for higher values of p (gate voltage) in the CDW 
regime. We observe that the boundary between CDW 
and RL2 phases as well as RL2-SC phase boundary shift 
slightly for higher values of V%. Fig. [3] shows the effect 
of NNN hopping on the phase diagram of Fig[2j We 
note that the RL2 phase and CDW region are squeezed 
on introducing the NNN hopping while expanding the 
region of SC phase. In Figs. [2]and[3l we have presented 
the explicit parameter dependence of QPD which is not 
possible to obtain from the Abelian bosonization study. 

If we consider only the interaction terms of the Hamil- 
tonian H%, for both Ai and A2 limits we get the con- 
dition Ezi < — ^ 1 for the boundary between RL1 and 
CDW. This condition is unphysical; we know experimen- 
tally that Eji, Ej2, Ezi and Ezi ar e all positive 
So the interaction part of Hamiltonian Hi alone is not 
sufficient to produce the whole phase diagram of Fig. 
[1] It indicates that an extended interaction term is re- 
quired to get the correct phase diagram. We focus phase 
boundary between RL2 and SC phase of Hi for A = Ai 
and A2, similar analysis of LL parameter gives the con- 
dition Ezi < which is again unphysical. This phase 
boundary also requires NNN hopping and NNN coulomb 
interactions. 

Here we do the analysis of the Hamiltonian H2 , which 
extends the range of hopping as well interactions of 
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FIG. 5: ti vs. V2 phase diagram for (a) p — 1.0 and (b) 
p — 0.5 (BKT) transition for two values of Vi = 0.4 and 0.3. 
Different points mark the different BKT points. Inset in (a) 
shows effect of finite values of t%. 



Hi, to NNN. Analysis for RL1 and CDW boundary, 
yields a quadratic equation in Ej%. The real roots of 
the quadratic equation do not impose negative sign for 
Ej\,Ej2, Ezi and Ezi unlike in the case with H%. 
Therefore, we can find physical values of these param- 
eters for which RL1 and CDW phase boundary can ex- 
ist. Similar analysis for the phase boundary between RL2 
and SC phase indicates that H2 can also support such a 
phase boundary. E>2 phase of Fig. [T]can be obtained only 
from the analysis of Hamiltonian H2 when E Z 2 exceeds 
some critical value like in the Majumdar-Ghosh model 

The above analysis is new and a consequence of 
extending the range of hopping and interactions. In 
Fig. [TJ we find two multi-critical point, P±, and P2. 
Mott, RL1, RL2, and superconducting phases coincide 
at the point Pi; RL1, RL2, and B\ or E>2 coincide at 
the P2 point. In our QPD, we obtained multi-critical 
points not only for the presence of different kinds of 
interactions in the SQD but also on application of 
external gate voltage on the SQD. In our QPD, the class 
of transition is changing as one moves away from the 
tip of the lobe in the Fig. [1] Our analysis allows to 
extend the classification of two types of Mott transition 
in the context of SQD. At the tip of each lobe, the 
system at fixed density is driven from the SC phase to 
insulating phase by changing the interaction parameter 
and this is Mott-U (Berezinskii-Kosterlitz-Thouless) 
transition. Away from the tip of the lobe, a change in 
the gate voltage drives the system from a SC phase to 
an insulating phase and this is the Mott-<5 transition. 



Fig. |4] shows that density-density correlation function 
F(j) = (< n m rij > — < Tij >< n m >) from the middle 
site of the chain, m. We study F(j) for different values 
of V2 and t 2 - We observe that as V2 increases the system 
crosses over from the 2K F CDW state to AK F CDW 
state. Increasing leads to decrease in CDW amplitude. 
The prediction of 4Kp CDW state is entirely new for this 
Hamiltonian. 

In Fig. O we present the DMRG analysis of BKT 
transition for p = 0.5 and p = 1.0 for different values of 
V\. We observe that for p — 0.5, SC phase occurs for 
small values of t\ compared to p = 1.0. For p = 0.5, 
CDW phase and RL2 boundary appears at t\ = 0.11 ± 
0.03 for Vi = 0.3 and 0.4 for different values of V 2 . We 
also note that £2 = 0.1 shifts the MI-SC boundary to a 
much smaller value of t\(< 0.03). 

In summery, we have presented the phase diagram of 
EBH model. Description of repulsive luttinger liquid 
phase regions and the AKp CDW states for higher values 
of V2 is entirely new. 
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